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TECHNICAL MEMORANDUM 


PERFORMANCE OF PHOTOMULTIPLIER TUBES AND SODIUM 
IODIDE SCINTILLATION DETECTOR SYSTEMS 

I. INTRODUCTION 


This report presents a discussion of the performance of photomulti- 
plier tubes (PMT's) and Nal scintillation counters. The test procedures 
and described results were performed to evaluate detectors for use on 
the Burst and Transient Source Experiment (BATSE) on the Gamma Ray 
Observatory (GRO). These detectors consist of 50.8 cm (20 in.) diameter 
by 1.27 cm (0.5 in.) thick Nal crystals viewed by three photomultiplier 
tubes. An important goal of the tests was to determine the optimum shape 
for a light collector housing to optically couple the PMT’s to the crystal. 

Section II presents a general discussion of performance characteris- 
tics of photomultiplier tubes, Nal scintillators, and light collector housings. 
This section is intended to be of general interest to any user of scintilla- 
tion detector systems. The following sections are more specific to our 
particular application. Section 111 describes procedures for measuring 
various parameters of detector systems, such as sensitivity, linearity, 
and energy resolution. Section IV presents the results of the tests for 
BATSE detectors. A brief summary of the conclusions is presented in 
Section V. 

This report is concerned primarily with linearity and variance of 
detector pulse height distributions. Topics not covered include speed of 
PMT’s, effects of magnetic fields, and temperature effects. 


II. PERFORMANCE CHARACTERISTICS 


An important parameter in these investigations is the relative vari- 
ance V, which is a measure of the energy resolution of the detector. For 
any distribution in x, the relative variance is defined by 



where x is the average value of x and a is the standard deviation (root 

mean square deviation) of the distribution. The advantages in using 
relative variance to describe the width of a distribution are that it is 




dimensionless and that it is additive when independent sources of varia- 
tion are combined. Scintillation crystal manufacturers usually specify 
energy resolution R in terms of the full width at half maximum (FWHM) 
of a spectral line, 


FWHM. 


R = 


generally expressed in percent. Here x is energy or PH A channel 
number. If the distribution is Gaussian, the variance can be obtained 
from the FWHM using the relation 


FWHM^ 
°x "" 2.354 


Thus, 


V = 0.1805 R^ 


A. Photomultiplier Tubes 

1. Variance . The variance of the pulse height distribution pro- 
duced by a photomultiplier tube illuminated by constant amplitude light 
pulses is a result of statistical fluctuations in the number of photoelectrons 
produced and statistical fluctuations in the multiplication by the dynodes. 
Assume that a light pulse deposits N photons on the photocathode of a 
PMT. A fraction the quantum efficiency, of these photons will eject 

photoelectrons. A fraction of these photoelectrons will be collected 

at the first dynode. Then n = ^ is number of photoelectrons 

collected at the first dynode. Since and represent Bernoulli trials, 

n will have a Poison jiistjpibution if N does . The variance will then be 

V = 1/n, where n = Nc^c . (the bar represents average values). The 
n 9 d 

variance in the output pulses, called the statistical variance V , is given 

by ® 


1 + V 


V = 
s 


m 


n 


2 


where the multiplication variance V accounts for statistical fluctuations 

m 

in the multiplication at the dynodes. V has been calculated for a few 

m 

idealized cases. If the multiplication at the dynodes is described as a 
Poisson distribution (an assumption found to be overly optimistic for most 
PMT's), then 



where 6 is the multiplication at the first dynode. If single photoelectron 

events produce an exponential pulse height distribution, then V =1. 

m 

Measurements of have ranged from 0.4 to 1 (Breitenberger [1] Swank 
and Buck ( 2] , Prescott and Takhar [ 3] ) . 

2. Gain . At normal operating voltages, the gain 5 of a single 

• n 7 

dynode varies with interdynode voltage approximately as V . In a 

PMT with k stages, with equal voltage per stage, the tube gain varies 

0 7k 

approximately as V ’ . Since PMT’s may have up to 13 stages, the 

output is very sensitive to high voltage, and voltage stability is important 
to maintain good resolution, 

3. Linearity . At high anode currents, PMT output becomes non- 
linear due to space-charge effects, usually strongest between the last 
two dynodes. An idealized derivation predicts that the limiting current 
density I between two electrodes is 


2.34 V 


3/2 


U A /cm' 


where d is the distance in centimeters between the electrodes and V is 
the potential difference in volts (Parsons [4]). Space-charge effects 
generally become significant at approximately 1 to 10 niA anode current. 
This value refers to pulsed current and is typically significantly larger 
than the allowable continuous current. Space charge cflects can be 
minimized by distributing the high voltage so that a relatively large 
potential difference appears across the last two dynodes. 

Nonlinearity may also bo present at high count rates if the dynode 
voltages are supplied by a resistor bleeder string. As current through 
the PMT increases, the voltages at the dynodes change. The redistribu- 
tion of voltage results at first in an increase in gain. At very high rates, 
the required current cannot be applied and the gain drops. This non- 
linearity has been examined analytically by Lush (51 and Land [6]. As 
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a rule of thumb, resistor bleeder string current should be at least ten 
times the average anode current. Capacitors placed across the last 
several dynodes can maintain linearity during individual large pulses or 
short bursts of pulses. Other techniques include using Zener diodes or 
transistors in the bleeder string and active control of the high voltage 
(Kovash and Blatt [7], Kerns [8], Hiebert, et al. [9]). 

Other causes of nonlinearity are cathode resistivity and fatigue 
at high light intensities (Keene [10]). These effects are negligible for 
most scintillation detector applications. 


B . Sodium Iodide Scintillators 

1. Luminous Efficiency . The luminous efficiency of a scintillator 
may be expressed as the ratio of light energy produced to ionization 
energy deposited or as the number of optical photons produced per unit 
ionization energy loss. Estimates of Nal luminous efficiency have ranged 
from 20 photons/keV to 50 photons/keV (Lindow, et al. [11], Bellian and 
Dayton [ 12] ; Akimov [ 13] , Knoll [ 14] ) . 

2. Linearity . It has been found that scintillation light output 

is not proportional to ionization energy loss in sodium iodide scintillators. 
Figure 1 shows light intensity per unit energy deposited as a function 
of gamma-ray energy (Zerby, et al. [15]). The curve was determined 
by a Monte-Carlo calculation. The points represent their measurements 
and can be seen to fit the predictions well. The output is approximately 
15 percent high at energies between 10 keV and 100 keV, relative to 
output at 662 keV. 

3. Intrinsic Line Width. The nonlinearity of sodium iodide intro- 
duces an additional source of variance. A primary gamma ray can produce 
secondary electrons in Nal in any number of ways, depending on the 
number of Compton scatters and the energy deposited per scatter. Each 
distribution of secondary electrons results in a different average light 
output due to the nonlinearity of the Nal. Thus, the distribution of 
light output for a monoenergetic beam of gamma rays is wider than a 
Poisson distribution of photons. The intrinsic line width V. is defined 

by ^ 


where is the total variance in the distribution of the number x of 

photons produced. For a Poisson distribution, V = 1/x and Vj = 0. 

Thus, ordinary statistical fluctuations in x are not included in the defini- 
tion of intrinsic line width. Vj has been calculated by Iredale [16], 

Zerby, et al. [15], and Narayan and Prescott [17]. The Monte-Carlo 
4 






E^. GAMMA-RAY ENERGY (ktV) 

Figure 1. Nonlinearity of Nal (Tt) output as a function of 
gamma -ray energy', from Zerby et al. [15]. The 
vertical axis is normalized to 1 at 662 keV. 

calculations of Narayan and Prescott [17] are shown in Figure 2, where 
Vj is plotted as a function of the energy of the incident photon. Vj is 

small at low energies, where photoelectric absorption dominates, and lises 
to a maximum at approximately SCO keV, where one or more Compton 
scatters are likely. At still higher energies, Vj drops due to the averag- 
ing effect of many Compton scatters. The step at 33 keV is due to the 
K edge of iodine. 


C. Transfer Variance and Light Collector Housings 

The probability that an optical photon emitted by the scintillator 
will result in a photoelectror jollected at the first dynode of a PMT is 
called the transfer efficiency T. T can be expressed as 


= 'r/q'd ■ 
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\ 



Figure 2. Intrinsic lino width of Nal (TO due to nonlinearity, 
as calculated by Narayan and Prescott ll7j. 

where is the efficiency for collecting light at the photocathode, 
is the quantum efficiency of the photocathode, and is the efficiency 
for collecting photoelectrons at the first dynode. 

The transfer efficiency will probably depend on the location wiihin 
the crystal at which a photon is emitted. The variance in the distribu- 
tion of transfer efficiencies is called the transfer variance V,j,. 

Large scintillators require some type of light collector housing or 
light pipes to bring the light to the photomultiplier tube. The design 
of the light collectors attempts to satisfy two ro(j ui remen ts: 

1) The efficiency of light collection must be high to keep the 
statistical variance low . 
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2) The uniformity of light collection must be high to keep the 
transfer variance low. 

These two goals can conflict in the geometrical design of a light collector 
housing. 

2 

Clay and Gregory [18] have studied light collection from aim 
plastic scintillator onto a single photomultiplier tube. They found that 
a number of white paints are satisfactory for the reflecting surfaces 
inside the housing and that uniformity improves with increasing distance 
between scintillator and PMT. The uniformity can be improved by placing 
a reflecting disc midway between the scintillator and the PMT to block 
some of the light that would fall directly on the PMT, but this procedure 
reduces the pulse height scxnewhat. 


D. Resolution of a Detector System 

The various sources of variance in a sodium iodide detector system 
can be combined to yield an output pulse variance for monochromatic 
incident gamma rays [1], 


Vo = . (1 . V,,) 



rn 

X T 


where 


X = average number of photons produced 

- variance in number of photons produced 

T - overage transfer efficiency 

variance in T (transfer variance) 

V = variance in multiplication of PMT 
m 

xT - number of photoelectrons collected at first dynode. 

Since is usually much less than 1 . the formula .nay be approximated 
as 
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where V„ is the transfer variance, V. = V - i is the intrinsic line 

_ _ X 

width, and = (1 + Vjj^)/(x T) is the statistical variance. 

III. PROCEDURES 


This section describes the test procedures in sufficient detail to 
enable the render to duplicate the tests. This section alsj provides the 
information necessary to judge the validity and accuracy of the results 
presented in the next section. 


A . Photomultiplier Tube Tests 

The photomultiplier tubes tested were EMI type D302B . which a; e 
11-stage tubes with 12.7 cm (5 in.) diameter bialkali photoenthodos. Also 
tested was an EMI type 9709R which has nine stages and an Sll photo- 
cathode. The dynodc resistor bleeder string used is shown in Figures. 
The resistances are those recommended by EMI for very high pulsed cur- 
rent NaKTO application. A relatively large resistance is used between 
the last two dynodos to reduce space-charge limiting of large pulses. 

EMI recommends that the resistor between the cathode and first dynode 
be selected to provide a potential difference of 450 V. We used 3.9 Mf^ 
for this resistance and tested other values. The tubes were enclo.sed in 
magnetic shields for all of the tests. 

Photomultiplier outputs were fed into an Ortec Model 113 preamplifier 
set at 0 pF input capacitance. Preamp pulses were fed into a TC 203HLR 
spectroscopy amplif ;r with > - 1 ns. BLR - out. The monopolar ac 
inverted output was fed into a Tracer Northcri! 1710 pulse height analyzer. 
The fine gain of the amplifier was set at 1.0, allowing direct comparison 
of pulse heights obtained at different coarse gains. A block diagram of 
the experimental setup is shown in Figure 4. 

PMT performance was tested using a green light-emitting diode (LED) 
covered by a translucent diffuser. Pulses were supplied by a BNC Bll 1 
tail pulse generator. The rise time was 0,5 us and the fall time was 
1.0 uS. The tail pulse generator was triggered by the 1 kc calibration 
pulse from a Tektronix 454 o.scilloscope to ensure reproducibility. The 
tail pulse generator was also used to gate the PHA, thereby greatly 
reducing noise and allowing the lower level discriminator on the PHA to 
be set to zero. The distance between the LED and the photo<'athode is 
not critical but should bo large enough to ensure fairly uniform illumina- 
tion of the photocathode (several inches). The I.ED electronic setup is 
also shown in Figure 4. 
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Figure 3. Resistor bleeder string used on 11-stage EMI 
Type D302B photomultiplier tubes. 


It is important to recognize that LED light output cannot be expected 
to be independent of pulse rate, since output is temperature dependent, 
nor is light output proportional to input pulse voltage. However, the 
number of photons emitted does appear to be Poisson distributed for con- 
stant input pulse heights. 
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Figure 4. Setup for testing photomultiplier tubes with 

LED pulses. 

1. Single Pliotoelectron Response. The LED can be used to deter- 
mine the response of the PMT to single photoelectrons. From this infor- 
mation, pulse heights can be converted to number of photoelectrons and 

can be measured. The procedure, in general outline, is as follows: 

1) Determine noise spectrum with LED off. 

2) Turn LED to an amplitude such that pulses are occasionally seen. 

3) Subtract noise spectrum, corrected for deadtime. 

4) Compute x^, the average channel number for single photoelec- 
trons, correcting for events in which more than one photoelectron was 
obtained . 

A detailed description of the procedure is given in the Appendix. 

2. Multiplieation Variance. Once a conversion between pulse height 
and number of photoelectrons is obtained, the PMT multiplication variance 

can be determined by measuring the variance of the pulse height dis- 
tribution as a function of number of photoelectrons. Details of the pro- 
cedure are given in the Appendix. should be measured at several 

different LED intensities. At high intensities, space -charge limiting may 
reveal itself as an apparent decrease in . 

The measurement of x^ need be made at only one value of high 
voltage. To determine x^ as a function of voltage, simply find relative 
PMT gain as a function of high voltage using a constant amplitude LED 
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pulse. To ensure that space-charge limiting does not occur at the higher 
voltages, verify that the relative variance at fixed LED amplitude does 
not drop with increasing voltage. At low voltages ( approximately 600 V) 
the variance increases sharp 'ly due to poor collection of the photoelectrons. 

Note that even if x, and V_ are not accurately determined, the 

1 m 

measured variance is still a valid measurement of the statistical variance 

V versus pulse height; When the scintillation crystals are tested, the 
s 

contribution of V to the total relative variance is, therefore, known quite 
accurately . 

3. Space-Charge Limiting . Nonlinearity due to space-charge effects 
was measured using several tecnniques. Several difficulties are involved 
in making such measurements. One difficulty is that we are interested in 
nonlinearity of pulses, not the dc current levels often used in investigations 
of nonlinearity. With pulses from sodium iodide, the PMT anode current 
rises rapidly to a peak value and then decays exponentially with a 
0.25 us time constant. The peak of the pulse will exhibit saturation before 
the exponential tail, as can be seen easily with an oscilloscope. Conse- 
quently, the pulse height nonlinearity depends on the timing characteris- 
tics of the output circuitry and the integration time of the amplifier. A 
fast peak-and-hold technique will exhibit greater nonlinearity than integra- 
tion of the pulse over several time constants. Another difficulty is spec- 
ification of the independent variable to avoid confusion. Traditionally, 
anode current is the independent variable, but this is precisely the quan- 
tity that saturates. We use "expected peak anode current", defined as 
the anode current, at the peak of the pulse, that would be expected in 
the absence of saturation effects. 

The advantage of using anode current as the independent variable 
is that the nonlinearity curves will not vary much from tube to tube nor 
show a high dependence on tube gain , although there is a dependence 
on voltage at the saturating stage. The anode current is measured by 
feeding the output of the PMT to a 50 '2 termination instead of the preamp. 
Since the RC time constant of the output circuit is so short (a few nano- 
seconds) , the voltage across the terminator tracks the anode current 
reasonably accurately. Therefore, the peak anode current equals the 
peak voltage across the terminator, measured from an oscilloscope, divided 
by termination resistance. 

Three different procedures were used to measure PMT nonlinearity. 

In the first, LED and a single optical filter were used. The LED was 
pulsed using a tail pulse generator with the rise and fall times set to 
match the pulse shape from Nal as closely as possible. The transmission 
of the filter was obtained by measurements at low LED intensities. The 
PMT output was determined with and without the filter for a wide range 
of LED intensities. Figure 5 shows the e.xperimcntal configuration. A 
difficulty with this procedure is that at high light intensities the anode 
current shows saturation effects even with the filter in place. Therefore, 
the linearity curve must be built up a step at a time, starting at low 
light intensities. It should be noted that LED nonlinearity precludes 
meaningful comparisons of pulse heights obtained at different input pulse 
amplitudes. 
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Figure 5. Setup for testing photomultiplier linearity at 
large pulse amplitudes. 


The second method is similar to the first except that the LED is 
set at its highest intensity and a series ot Wratten filters ranging from 
ND 0.5 to ND 4 are interposed. This procedure is faster than the first 
but requires accurate knowledge of the filter densities. Cax*e must be 
taken that no light can leak around the filter to the PMT. 

The third method employs radioactive sources and Nal crystals to 
produce the light pulses. The ratio of the pulse heights produced by 
two radioactive sources is determined as a function of PMT high voltage. 
The ratio of the light input to the PMT is not ; cessarily the ratio of 
the gamma-ray energies, since Nal output is not linear with energy 
deposited. The input ratio is determined by measurements at low values 
of PMT high voltage. With this procedure, the linearity curve does not 
apply to any single value of high voltage. However, the gain of an 
11-stage tube is so strongly dependent on voltage that the entire linearity 
curve can be produced with a variation of approximately 300 V. Advan- 
tages of this toclinique arc that there is no question about the pulse 
shape accurately matching that of sodium iodide and that the nonlinearity 
can be determined at particular values of gamma- ray energy for the 
specific detector system being tested. 

High-Rate Effec ts. The gain of a PMT using a resistive dynode 
bleeder string depends on count rate, as described in Section II. A. This 
effect cannot be tested reliably with an LED, since LED light intensity is 
not necessarily independent of pulse rate. Instead, the PMT was coupled 
to a 12.7 by 1.27 cm (5 by 0.5 in.) Nal crystal illuminated by 60 keV 

241 

gamma rays from a 10 mCi Am source. The PMT output was fed to a 
50 termination resistor. Average anode current was determined from 

the equation 
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where V is the peak voltage across the 50 Q resistor when a 60 keV gamma 
ray is detected, t = 0.24 ys is the Nal scintillation decay time, r is the 
count rate, and R = 50 This equation may be understood in the follow- 
ing way: V /R is the peak anode current; Vt/R is the integrated current 

during a pulse, or the charge transferred in one pulse; I^, the average 

rate of charge flow, is then (Vx/R) r. The rate r is varied by varying 

241 

the distance between the Am source and the crystal. The 60 keV 
pulse height is then measured as a function of count rate. The average 
anode current is measured only once, at low count rate. The expected 
average anode current is then scaled as the count rate. The actual 
average anode current is just the expected average anode current times 
the fractional change in pulse height. It is useful to compare average 
anode current to the bleeder string current at zero count rate. With high 
voltage set at 1600 V, the quiescent bleeder string current is 80 pA. 

5. Other Tests. PMT gain and resolution were determined as a 
function of high voltage using LED pulses. Resolution was investigated 
as a function of cathode to first dynode voltage by changing the appro- 
priate resistor in the bleeder string. PMT noise rate was measured as 
a function of time since exposure to room lights. 


B. 12.7 by 1.27 cm (5 by 0.5 in.) NaI(T!L) Crystal 

After determining the conversion between pulse height and number 
of photoelectrons and determining for a photomultiplier tub?, addi- 
tional sources of variance in a NaI(T£) detector system were investigated. 
A Harshaw 12.7 by 1.27 cm (5 by 0.5 in.) NaI(T£) crystal with a reflect- 
ing rear surface was coupled to a 12.7 cm (5 in.) photomultiplier tube 
using Dow Corning optical coupling compound. Gamma-ray sources of 
known energy were used to determine pulse height (and thus n ) and 

total variance as a function of energy deposited in the crystal. Tests 
were made with uniform illumination of the crystal by the gamma rays 
and with the gamma rays collimated by a lead absorber to a small region 
in the center of the crystal. The difference in resolution between colli- 
mated and uncollimated beams provides an estimate of the transfer variance 
in the uncollimated case. 

The gamma-ray sources used and the energies emitted are given in 
57 

Table 1. Co was not used because this isotope emits gamma-ray lines 
at 122 keV and 136 keV, resulting in a single anomalously broad line at 
approximately 124 keV. Many detector manufacturers specify resolution 

137 

at the 662 keV line of Cs . As will be shown in Section IV, however, 
the variance at 662 keV is dominated by Nal intrinsic line width. There- 
fore, resolution at this energy does not provide a sensitive indicator of 
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TABLE 1. GAMMA-RAY SOURCES 


Source 

Energy (keV) 

Mo K X-ray 

17.4 


Ag K X-ray 

22.1 




Variable Energy X-ray Source 

Ba K X-ray 

32.2 


Tb K X-ray 

44.5 



59.5 


22.1 


88 


30.6 


81 


356 

Na22 

511 


1275 

Cs^'^ 

662 

Mn^' 

835 


light collection efficiency or uniformity. In using sources emitting gamma 
rays of more than one energy, the lower energy peaks may be distorted 
due to background from lines of higher energy. Background spectra 
were subtracted from source spectra before deter ..ining the peak and 
width of gamma-ray lines. 


C. Light Collector Housings and 
50.8 by 1.27 cm (20 by 0.5 in.) NaI(TO Scintillation Crystals 

A number of tests were made to determine the characteristics of 
detectors using 50.8 by 1.27 cm (20 by 0.5 in.) crystals viewed by three 
PMT's. The main goals of these tests were to determine the light emission 
efficiency and uniformity of several types of crystals and to determine 
the optimum geometry for the light collection housing. Figure 6 illustrates 
the geometry of the detector system. The crystals are viewed by three 
photomultiplier tubes on a circular plate. A cone or cylinder supports 
the PMT plate above the crystal. The interior surfaces of the cone and 
PlvlT plate were painted with Eastman No. 6080 white reflective coating, 
applied approximately 1 mm thick in accordance with the manufacturer's 
recommended procedures. 
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Figure 6. Light collection housing mechanical configuration for 
testing 50.8 by 1.27 cm (20 by 0.5 in.) NaI(T?,) crystals. 

The anodes of the three PMT's were summed using a resistance 
network. Separate high-voltage power supplies were provided for each 
tube and were adjusted to equalize the gains. Figure 7 illustrates the 
experimental arrangement. The LED was used to balance tube gains and 
to determine the single photoelectron pulse height and multiplication vari- 
ance for the summed output using the same procedures as for individual 
tubes. Collimated and uncollimated gamma- ray beams were used to investi- 
gate transfer variance. 



Figure 7. Electronic setup for testing 50.8 by 1.27 cm 
(20 by 0.5 in.) crystals using three PMT's. 



IV. RESULTS 


A. Photomultiplier Tube 

This section describes the results of tests performed on EMI type 
D302B photomultiplier tubes. The procedures for some of the tests are 
described in Section III. A, A few of the tests were performed on three 
different tubes. 

1. Noise. The PMT count rate as a function of time after exposure 
to room lights was investigated. Figure 8 shows a plot of count rate ver- 
sus time. The rate decreases rapidly in the first few hours. Subse- 
quently, there is a slow decrease to a minimum, followed by a slow 
increase. This slow variation is probably due to the change in ambient 
temperature throughout the night. This test shows that photomultiplier 
tubes must be kept in darkness at least several hours before performing 
any tests in which small output puLses must be detected, as in determining 
the single photoelectron response. 

2. Cathode to First Dynode Voltage . EMI recommends that the 
voltage between the cathode and the first dynode of its 12.7 cm (5 Li.) 
tubes be 450 V. Resolution of an LED pulse of fixed intensity was meas- 
ured for values of cathode to first dynode voltage V„ ranging from 

200 to 508 V. The results are shown in Table 2, where V . is the 

K “A 

cathode to anode voltage. The LED intensity was such that each pulse 
yielded approximately 55 photoelectrons. The resolution is independent 
of jjj, indicating that this voltage is not critical within the range 

examined . 

3. Single Photoclectron Spectrum. As described in Section III.B, 
an approximation to the single photoelectron spect»'um can be obtained by 
pulsing the LED at an amplitude low enough that the average number of 
photoelectrons per pulse is much less than one. Using this procedure, 
the single photoelectron spectrum for EMI tubes was found. A typical 
spectrum, shovvn schematically in Figure 9, is approximately exponential 
with a flattening at the lowest pulse heights. 

Multipl icatio n V^a ria nce . Multiplication variance V^ is determined 

by measuring the variance as a function of number of photoelectrons, as 
described in Section lll.B. Since V - ( 1 + V' ) /n , it is customary to 

plot variance versus (1/n^^). A straight line of slope (1 + V^) is 

expected. Figure 10a shows such a plot. A value V^ = 0.64 fits the 

data well for the three PMT's tested. No variation of V^^ with voltage is 

evident. Figure 10b provides an expanded view of the region of Figure 
10a near the origin. 


RATE (COUNTS/HR) 



Figure 8. Photomultiplier noise rate versus time since exposure to 

room lights. 


TABLE 2. CATHODE TO FIRST DYNODE VOLTAGE 


V 

K-Dl 

(Volts) 

^K-A 

(Volts) 

Resolution 

(FWHM) 

200 

1000 

40^ 

263 

1000 

40% 

339 

1000 

39% 

508 

1600 

40% 
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Figure 9. Idealized shape of single photoolectron spectrum. 

5. Gain versus High Voltag e. Photomultiplier gain as a function 
of high voltage is plotted on a log-log scale in Figure 11 for an 11- stage 
and a 9- stage tube. The data are well represented by a .straight line. 

The slope is 7.6 for tne 11- stage tube and 5.5 for the 9- stage tube. 

These values for the slope imply a gain per stage proportional to (HV)^, 
where \ = 0.6 to 0.7, about as expected. At higher voltages, the gains 
decrease due to space-charge limiting. 

6. Space-Charge Limiting . Nonlinearity due to space-charge limit- 
ing was measured by the methods described in Section III. A. 3. The 
results are presented in Figure 12. Nonlinearity, expressed as a percent- 
age, is defined as follows: 

Expected Pulse Height - Measured Pulse Height 

Nonlinennty = Eipgcro ? P^ lso lloiRlH ' 
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a. Relative variance of PMT output versus 1/n 

pe 

multiplication variance is approximately 0.64. 


The 


Figure 10. Relative variance. 
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UM VOLTAOl 

Figure 11. Log pulse height versus log high voltnge for 11-stnge 

and 9- stage PMT's. 






The expected peak anode current is the peak anode current that would 
be obtained without space-charge effects. The different symbols repre- 
sent data taken using various combinations of PMT, high voltage, and 
technique. The dashed line labeled EMI represents the expected perform- 
ance presented in the EMI Photomultipliers catalog. 

The data indicate that the various techniques agree reasonably well 
and that the nonlinearity is approximately 5 percent at 1 mA peak anode 
current and approximately 25 percent at 10 mA peak anode current. The 
PMT's performed significantly worse than the EMI predictions, which 
provided no details on the measurement technique. 

7. Nonlinearity at High Kates . Nonlinearity occurs at high count- 
ing rates as a result of large average current flow causing redistribution 
of the dynode voltages. Figure 13 shows a plot of measured average anode 
current I as a function of average light intensity, which was varied by 

• ** 241 

changing the distance of an Am source from the 12.7 by 1.27 cm (5 
by 0.5 in.) crystal. The quiescent current through the resistor bleeder 
string Ij^ was 80 y A. At low anode currents the measured tube nonlinear- 
ity can be well represented by AGain/Gain = 0.6 (I /ly.) , in reasonable 

d o 

agreement with the calculations of Land [6] and Lush [5]. Therefore, 
to maintain linearity within 1 percent, the bleeder string current must 
be 60 times the average anode current. The tube saturates at an anode 
current of approximately 0.55 Ij^. 

8. Resolution versus High Voltage. The resolution of the D302B 
photomultiplier tubes was examined to see if it depended on high voltage. 
LED pulses yielding about 65 photoelectrons were used. The results are 
shown in Figure 14. The resolution is approximately constant at approx- 
imately 37 percent over the range 700 to 1500 V. At lower voltages, the 
resolution rapidly degrades, presumably due to poor focusing of photo- 
electrons. At higher voltages, the resolution appears to improve slightly. 
This may be due in part to a decrease in V and in part to space-charge 
limiting. 


B . Sodium Iodide 

1. Luminous Efficiency . Although the luminous efficiency of Nal 
was not measured directly, a more useful parameter, the number of photo- 
electrons collected at the first dynode per keV of energy' deposited in the 
crystal, was determined. We call this parameter the system efficiency e. 

It may be expressed as c = (X T)/E, where X is the average number of 
photons produced by absorption of a gamma ray of energy E , and T is 
the average transfer variance. Since c depends on transfer efficiency 
and gamma-ray energy, the following conditions for its measurement are 
specified : 

1) The energy is 662 keV (Cs^^^ line). 

2) The crystal is a Harshaw 12.7 by 1.27 cm (5 by 0.5 in.), 
directly coupled to the PMT using Dow Corning optical coupling compound. 
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LIGHT INTENSITY (ARBITRARY UNITS) 

Figure 13. Deviation from linearity due to high rates. is the 

quiescent current through the bleeder string and equals 

80 It A. The dashed line represents the relationship 

A Gain /Gain = 0.6 I /I, , where I is the 
a b a 

anode current. 

3) The gamma- ray beam is collimated to a spot at the center of 
the crystal. 

Under these conditions, r was determined to be 13.9 photoelectrons /keV. 


A lower limit to the Nal luminous efficiency is found by assuming 
perfect light collection = 1), perfect photoelectron collection (f'^j 1) . 

and peak quantum efficiency ^ = 0.27 for a bialkali photocathode. The 

luminous efficiency is > 50 photons/keV, which agrees well with the 
Harshaw specification and corresponds to 15 percent conversion efficiency 
at 662 keV. 


The parameter f, depends upon the measurement of the pul.se height 
for single photoelectrons and thus propagates errors in that measurement. 
However, a lower limit to c can be obtained by measuring the resolution 
at a low energy and attributing all of the variance to photoelectron 
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Figure 14. Resolution versus high voltage. The large increase below 
700 V is probably due to poor collection of photoelectrons. The 
decrease at voltage greater than 1500 is probably caused 
by a combination of space -charge limiting and a 

decrease in V . 

m 

statistics (i.e., neglect transfer variance and multiplication variance). 

941 

The resolution at 60 keV (Am ) is 11.8 percent FWHM, indicating 
e > 6.6 photoelectrons /keV , and luminous efficiency > 24 photons/keV. 

2. Nonlinearity. Measurements of pulse height versus gamma-ray 
energy were made to determine the nonlinearity of Nal(TC) output. The 
PMT high voltage was low (1000 V) to ensure that space-charge effects 
were negligible. Figure 15 shows a plot of light intensity per unit energy 
(normalized to 1 at 662 keV') versus gamma-ray energy. The curve is 
taken from Zerby, et al. [15]. The solid squares represent our observa- 
tions using a 12.7 by 1.27 cm (5 by 0.5 in.) Harshaw crystal directly 
coupled to an EMI 12.7 cm (5 in.) PMT using Dow Corning optical 
coupling compound. The beam of gamma rays was collimated to an approx- 
imately 1 cm spot at the center of the crystal face. Our observations are 
in excellent agreement with the observations and calculations of Zerby, 
et al, [15]. 


4 ^ 
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E^. GAMMA-RAY ENERGY (kaV) 

Figure 15. Measured nonlinearity of Nal(Tv) output versus gamma- 
ray energ>'. Results obtained here are shown as solid s(iuares 
and agree well with calculations and measurements of 
Zerby , et al. [ 151 . 

3. V'ariance. Variance of pulse height versus gamma -ray energy 
was also determTned during the measurements of nonlinearity described 
in the last section. The results are plotted in Figure 16. The solid 
circles represent data obtained with a collimated beam of gamma rays. 
The open circles represent data taken when the gamma sources were 
uncollimated and distant ciiough ( ' ] ni ) that the crystal may be consid 

ered to be uniformly illuminated. I'he solid line represents a calculation 
of variance from the etiuation 


1 + V 


V - V.^, + (I + 


where the intrinsic line width [V^ (l/x)| was taken from Narayan and 

Prescott [17], x T (- number of photoelectrons collected) and were 

m 

measured, and V,j, was sot ecpial to 0. Thus, the curve represents a 

minimum variance with no free parameters. The pt)ints for the collimated 
beam lie slightly above the calculated curve, imlicating an additional source 
of variance, presumably transfer variance. The discrepancy increases 
with decreasing energy. The uncollimated beam results in an additional 
transfer variance of approximately 0.002. The ratio of light collection 
efficiencies for uncollimated versus collimated beams is 0.95. 
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Filfure 16. Relative variance versus gamma-ray energy for collimated 
(solid circles) and uncollimated (open circles) gamma-ray beams 
incident on the 12.7 by 1.27 cm (5 by 0.5 in.) crystal. The 
line shows a calculation assuming transfer variance = 0. 


C. Light Collector Housings 

The procedure for testing light collector housings on the 50.8 cm 
(20 in.) crystals was described in Section III.C. In a preliminary test, 
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pulse height and variance of the 88 keV line of Cd ‘ were measured for 
six configurations of light collector using an off thc-shelf Bicron crystal. 
The six configurations arc shown In Figure 17 drawn to scale. Parameters 
specifying the geometry are the height h of the PMT plate above the 


L'V 


r 



scintillator, the radius r of the PMT plate, and the distance d between 
tube centers (Fig. 6). The results are presented in Table 3, where 
gfewnetry parameters, light collection efficiency, and transfer variance are 
tabulated. The light collection efficiency was determined by comparing 
the pulse height with that obtained using the 12.7 by 1.27 cm (5 by 
0.5 in.) crystal, assumed to have a light collection efficiency of 1. 
Uncollimated gamma-ray beams were used in all configurations, and col- 
limated beams were also used in most of the configurations. 


TABLE 3. LIGHT COLLECTOR HOUSING PERFORMANCE AT 

88 keV 


Confltrumtion 

Plato 

Height 

PUitc 

KitdiuM 

Tube 

Sepftrntioti 

I.IrM Cu’lection 
Kf(?oierw.'v T 

RrsdUituVi 

r Vurtimca 

Number 

h (cm) 

r ( cm 1 

d ( cm ) 

Cnllimatrit 

fT'iu-ollfe'iTiTr 




1 ncoHiniNted 

1 


2 \ 

2.1 

0.003 

O.ftOA 

0 40 

0.46 

0,017 

0.026 

2 

M 


J.5 


O.OKS 


0 71 


0.077 

3 

31 

10 

10 


0 000 

0 U 

0.44 

0.019 

0.021 

4 

17 

10 

10 

0 . n 

0. OHO 

0.44 

0 . b2 

0,024 

0.0.46 

5 

9 

10 

10 

«. u 

0 ll«.’ 

0 . 0.1 


0.060 

0. 107 

6 

17 


23 


6 0H8 


0 


0.0.18 


The light collection efficiencies are remarkably uniform from one 
configuration to the next, but very low compared to the 12.7 by 1.27 cm 
(5 by 0.5 in.) crystal. Because the different geometries had vastly dif- 
ferent reflecting surface area, the constancy of the light collection effi- 
ciencies implies that the Eastman No. 6080 high - reflectance coating had 
negligible absorption of light. The low efficiency of light collection 
(a- 9 percent) is apparently due to absorption in the crystal. [It is 
assumed that the luminous efficiency of the 50.8 by 1.27 cm (20 by 
0.5 in.) crystal is et|ujil to that of the 12.7 by 1.27 cm (5 by 0.5 in.) 
crystal]. In fact, a visual inspection of the crystJil supports this 
hypothesis. The 50.8 cm (20 in.) crystal is noticeably greyer than the 

12.7 cm (5 in.) crystals. Two further tests of the 50.8 cm crystal were 
performed. In one test, a PMT was directly coupled to the face of the 

50.8 cm crystal with Dow Corning optical coupling compound and illuminated 
with a collimated beam of 88 keV gamma rays. The light collection effi 
ciency was 0.53. .Although some of the light is expected to e.sc.ape in 

this geometry, the light collection efficiency is significantly less than 
that of the 12.7 cm crystal. In another lest, the 50.8 cm crystal was 
replaced by the 12.7 cm crystal placed at the center of a 50,8 cm disc 
of white styrofoam, using configuration 4 for the light collector hou.sing. 
The light collection efficiency increased from 0,11 to 0.31 for collimated 
gamma rays. 

The poor optical properties of the 50.8 cm (20 in.) crystals are 
probably due to absorption of light by the surf.jce of the aluminum con 
tainer. This surface has a sprayed coat of white paint, rather than the 
high- reflectance MgO surface of the 12.7 cm (5 in.) c:*;>' uiis. The optical 
properties of the large crystals are adequate for ti.cir intended applicf.tion 
as gamma- scan camera plates. 




The transfer variance is strongly dependent on light collector geom- 
etry. As Table 3 shows. ranges from 0.021 to 0.107. Figure 18 shows 

a plot of resolution versus distance h from the PMT's to the crystal for 
an uncollimated beam of 88 keV g^mma rays. The resolution improves with 
increasing distance, leveling off above approximately 20 cm. This depen- 
dence is due to a more uniform "view" of the crystal by photomultipliers 
at greater distance and implies that a significant fraction of the collected 
light is unreflected after emission from the crystal. 



Figure 18. Resolution of 88 keV line versus height of PMT 
top plate above crystal for the 50.8 by 1.27 cm 
(20 by 0.5 in.) crystal. 

After these tests, three additional 50.8 by 1.27 cm crystals wore 
obtained to study their optical properties in our application. The four 
crystals tested were; 
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A - Bicron off-the-shelf crystal with painted backing 

B - Bicron crystal with high -reflectance paper backing 

C - Harshaw crystal with high- reflectance paper backing 

D - Harshaw crystal with 0.25 in. Pyrex glass cemented to both 

sides. 

Crystal D was selected for study because it may provide greater mechani- 
cal strength. The window on one side of that crystal was covered with 
Millipore filter paper to act as a light-reflecting surface. 

The light outputs of the crystals were compared using a standard 
cylindrical light collection housing. The height h was 35 cm, plate radius 
r was 28 cm, and tube separation d was 30 cm. This geometry is not 
particularly efficient, but the relative light output of the crystals them- 
selves was the focus of this test. Table 4 presents the results using 
88 keV gamma rays. Crystal B exhibited the highest system efficiency 
and the best resolution, but all three of the new crystals were signifi- 
cantly better than the originally tested crystal (crystal A). 

Next, a series of tests were undertaken to determine the optimum 
light collection housing for the crystal with the best optical properties 
(crystal B). The results from the tests on crystal A are not directly 
applicable but provided guidance for choosing the designs that are most 
likely optimum. It should be noted that the best geometry represents a 
compromise between light collection efficiency and uniformity. Since the 
statistical variance decreases with energy, a specific design will be 
optimum at one energy only. We chose to design for best performance 
at 88 keV because the average energy of BATSE events should be near 
88 keV mid because a clean 88 keV gamma-ray peak is obtained from 
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Cd . The source was approximately 60 cm from the ci*ystal, giving 
uniform illumination to within approximately 20 percent. The configima- 
tions tested using crystal B are shown in Figure 19. drawn to scale. In 
configurations 1 and 2. the housing was a cylinder, permitting any height 
h, up to 36 cm, to be examined. 

Resolution was measured with all three tubes operating and with 
only two tubes operating. The latter test was made because there is a 
possibility of PMT or high voltage power supply failure during the 2 year 
GRO mission. The results are shown in Figure 20. where resolution is 
plotted against housing height h. The points below the dashed line at 
41 percent resolution represent data taken witli three PMT's operating; 
the points above the line are for two PMT's operating. The different 
symbols refer to tlie different configurations shown in Figui*e 19. The 
two curves show the trend of the data. 

The effect of tube separation d on re.solution can be soon by com- 
paring configuration 1 (d - 29 cm) with configuration 2 (d - 23 cm). 

With all three tubes working, the difference in resolution is slight; but 
with two tubes only, configuration 2 (with closer tubes) is superior 
(better resolution), especially at small heights, as might be expected 
intuitively. ... 


TABLE 4. COMPARISON OF OPTICAL PROPERTIES OF 50.8 cm 

(20 in.) CRYSTALS 


Crystal 

System Efficiency c 
(photoelectrons/keV) 

Resolution @ 88 keV 
(%) 

A 

0.76 

48 

B 

1.05 

38 

C 

0.91 

41 

D 

0.95 

41 


Figure 20 also shows that the resolution with two tubes degrades 
seriously for all configurations having height h less than approximately 
25 cm. The cone configurations provide better resolution than the cyl- 
inders at any value of h. The best configurations were numbers 3 and 
5, both tall cones. The best resolution was 32 percent for three tubes 
and 42 percent for two tubes. 

Based on these tests, the optimum geometry for the light collection 
housings is a cone, 25 to 30 cm high, and u tube separation of approxi- 
mately 20 cm. 


V. CONCLUSIONS 


The following is a brief summary of the findings of this study: 
EMI type D302B Photomultiplier Tubes 

• V - 0.64 

m 

• High voltage 700 V for good resolution 

• Cathode to first dvnode voltage not critical between 200 
and 500 V 

• Gain per stage > (HV)^ . where ■, = 0.6 to 0.7 

• Space-charge effects arc 5 percent at 1 mA peak anode 
current and 25 percent at 10 niA 

• Bleeder string current must bo 60 times average anode 
current for linearity w. bin 1 percent 


*■ 
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h (cm) 

CONFIGURATIOM NO 


1 - O 

7-0 

3 - 0 

4 - • 

6 - ■ 

Fipurc 20. Rc.solution versus height h of light collector housing 

using crystal B . 

ScMhum Iodide 

• With good light collection, 13.9 photoolectrons collected 
at first dynode per keV deposited (at 662 keV) 

• Nonlinearity in agreement with Zerbv et al. (151 




50.8 cm (20 in.) Crystals and Light Collection Housings 

• Bicron crystals with paper reflector had best light output 

• Best geometry is a cone of height 25-30 cm, with tubes 
^aced 'V 20 cm 

• With best geoineiry, energy resolution at 88 keV is 'v 32 
percent with three tubes operating, and 42 percent with 
cmly two tubes operating. 
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APPENDIX 


PROCEDURE FOR DETERMINING SINGLE PHOTOELECTRON 
SPECTRUM AND MULTIPLICATION VARIANCE 


The f(dlowing procedure is used to determine the spectrum of single 
photoeleotron events and to determine the average pulse height corre- 
sponciinfiV to one photoelectron: 

1) Set up electronics as shown in Figure 4. The PMT must be 
kept in darkness for several hours before making measurements. Set 
PHA lower level discriminator to zero and select coincidence mode. 

2) ih preamplifier on but with high voltage off, set amplifier 
gain as high as possible without allowing preamp noise to go above the 
first channe. or two of the PHA. Deadtime frwn preamp noise must be 
small. 


3) With LED amplitude at zero, turn on high voltage and increase 
until PMT noise (which is largely due to single photoelectron events) 
extends approximately 20 to 30 channels. With the EMI 11-stage tubes, 
approximately 1400 V was found to be satisfactory. 

4) Set the PHA livetime to some standard value (500 or 1000 s 
is satisfactory). Turn the LED amplitude sufficiently high to have an 
analyzed pulse every trigger. Acquire a spectrum. Integrate the spectrum 
to determine the total number of triggers per run N^. 

5) Acquire a spect:'um with the LED amplitude at zero to obtain 
the background spectrum. Determine the number of background events 
per run Ng and save the background spectrum. It is necessary that 

Ng $ 0.1 N^. This criterion may be difficult to achieve with PMT's having 

less than 11 stages. If the background rate is too high, try a higher 
voltage on the tube and a lower amplifier gain. If necessary, set lower 
level discriminator to cut out noise. High backgrounds may also be caused 
by taking measurements too soon after exposure of the PMT to room lights. 

6) While observing amplifier output on the oscilloscope, increase 
LED pulse amplitude until pulses are occasionally seen. However, keep 
amplitude low enough that most (> 80 percent) of the triggers show no 
pulse. This requirement ensures that most of the events observed are 
due to single photoelectrons. 

7) Acquire a spectrum , called the background + LED spectrum . 
Integrate to find the number of events analyzed N,. . 
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8) The spectrum of LED pulses is obtained by subtracting the 
background spectrum from the background + LED spectrum and correcting 
for deadtime. The background spectrum is corrected by multiplying by 
( 1 - fj^) , where fj^ is the fraction of triggers in which an LED pulse is 

present and is found from the equation 




B 

B 


This correction, which accounts for those triggers during which a back- 
ground and an LED pulse are present, can be made easily on the TN 1710 
PH A by using the multiply feature of the Data Processor Module. The 
validity of the background subtraction is assured if the first channel is 
high in the background spectrum and the background + LED spectrum due 
to preamp noise, but lies on a smooth curve in the LED spectrum. 

9) Compute the average channel number x of the LED spectrum. 

10) The LED spectrum is not equivalent to the single photoelectron 

spectrum, since there is occasionally more than one photoelectron per LED 

pulse. The average number of photoelectrons per analyzed pulse ri is 

pe 

determined from the known rate of triggers in which at least one photo- 
electron is present. From Poisson statistics, it can be shown that 


£n ( 1 - L ) 

n = - 2 = — 

pe r;; 

As fj^ ^ 0, 1, and the LED spectrum asymptotically approaches 

the single photoelectron spectrum. 

11) Determine the average channel for single photoelectrons 

from the equation x, - x In . Additional runs should be made at various 

LED pulse amplitudes to verify that is insensitive to amplitude. For 

f. << 1, statistical errors in x become large. As f ->■ 1, statistical errors 
Li L 

in n become large . 
pe ^ 

The multiplication variance can now be determined using the 
same experimental setup. The procedure is as follows: 
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1) Increase the pulse generator amplitude until an LED pulse is 
present every trigger (fj^ = 1). 

2) Acquire a spectrum and deierriine the variance V and mean 
channel number k of the pulse height distribution. For sufficiently large 
LED pulses, the spectrum is approximately Gaussian. In this case, the 
peak channel may be used for the mean, and the FWHM may be used to 
find the relative variance from the equation 


V = (FWHM)^ 

2.354 X 


3) From the mean channel number x determine the mean number of 
photoelectrons = x/x^. 


4) Determine the multiplication variance from 


V = V n 
m s pe 


1 
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